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I. INTRODUCTION 



Since the advent of gauge theories in the 1970-s there is a continous need for their 
application at high temperatures characteristic to high energy collisions of heavy nuclei or 
the early universe. Even asymptotically free theories, like QCD, prove to be far from trivial 
at high energy density or temperature. 

Rather recently a promising method has been invented, the hard thermal loop (HTL) 
resummation , for treating field theories at high temperature consistently. This method is 
analytical, complementing extended numerical simulations on a space time lattice at finite 
temperature. 

One of the most famous results of the HTL method, the gluon damping rate at zero mo- 
mentum ||, was found to coincide with half of the Lyapunov exponent of a classical real-time 
lattice description of nonabelian gluon fields ||. An attempt for understanding this at the 
first sight rather surprising agreement between two different quantities of high temperature 
gluon dynamics has been recently made by emphasizing that the main contribution to the 
gluon damping rate of order g 2 T stems from long wavelength gluons (i.e. from the infrared 
part of the gluon distribution) and hence can be a measure of the self-ergodizing speed of 
classically chaotic gluon fields ||]. 

In the present article we investigate whether a similar correspondence can be found in 
the coupled Yang-Mills Higgs-doublet system between the Higgs damping rate calculated 
in the HTL approximation and the Lyapunov exponent calculated on the lattice. If yes, it 
would let us surmise that something inherently classical is contained in the HTL method 
in the order of g 2 T and their results in principle could be simulated by classical dynamical 
calculations. If not, it would warn us that there is something uniquely particular for the pure 
gauge field systems. In fact chaos is opposite to damping, but because of time symmetry 
exponential growth and damping can have the same rate in given approximations. These 
are two facets of the same dynamics studied by using different boundary conditions in time. 

Although the Yang-Mills-Higgs system has already been studied classically and found to 
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be chaotic 0, for a throughout comparison with the new HTL results presented here we 
have to redo and extend that calculation, too. 

II. DAMPING RATES 

The aim of this section is the calculation of the Higgs and gauge boson damping rates at 
zero momentum and high temperature for comparing it with the Lyapunov exponent of the 
classical SU(2) Yang-Mills system coupled to a Higgs field ||. We start from the Lagrangian 

£ = ~\ FjL, F? + (D^V (D"$) + /i 2 $ f $ - A ($* $) 2 , (1) 

where $ is a charged Higgs doublet and = d^—g A%T a the covariant derivative containing 
the SU(2) generators T a . At high temperatures above the phase transition ~ T 2 and 
the quadratic term can be neglected besides the quartic term. Throughout this paper we 
consider fj, — 0. 

The damping rates follow from the imaginary part of the Higgs and gauge boson self 
energies, respectively. In order to obtain gauge invariant results, complete to leading order 
in the coupling constants, we adopt the Braaten-Pisarski method It is based on the 
distinction between hard momenta of the order of the temperature T and soft ones of the 
order gT or y/XT. The Braaten-Pisarski method amounts to the use of effective Green 
functions constructed by resumming the HTL diagrams. Calculating the gluon damping 
rate at zero momentum in this way solved the plasmon puzzle, a famous problem of finite 
temperature gauge theories @. 

The first step is to extract the HTL contributions to the effective Green functions, which 
arise from one-loop diagrams where the internal momenta are hard. Assuming all external 
momenta of the Green functions to be soft, the HTL corrections are of the same order in g 
and A as the bare Green functions. In the case of a Yang-Mills field coupled to the Higgs field 
described by the Lagrangian ([[J) there are no HTL corrections to effective vertices involving 
Higgs lines, as can be shown by power counting JTJ. This observation holds also for scalar 
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QED which is apart from the gauge boson self coupling and SU(2) factors identical to 
(0), as well as for the Yukawa theory ||. 

The HTL Higgs self energy, from which the effective Higgs propagator is constructed 
by resummation, follows from the diagrams of Fig.l, where the internal momenta of the 
polarization graph are hard. Adopting the usual Feynman rules ]7|], taking into account 
symmetry factors 1/2 for the tadpole diagrams, we obtain the HTL Higgs self energy using 
the Matsubara formalism, 

^ = <W (I A + T 2 (2) 

with the SU(2) indices a,/3 of the fundamental representation. As in the case of scalar QED 
H and the Yukawa theory |J the HTL self energy is momentum independent and real. 

The effective Higgs propagator is given by resumming the HTL self energy (Q) in a 
Dyson- Schwinger equation leading to 

where K 2 = k 2 — k 2 , k = |k| and m 2 H = (4A/3 + 3g 2 /16) T 2 is the square of the thermal 
Higgs mass. 

The HTL gauge boson self energy is shown in Fig.2. The first three diagrams, contain- 
ing the gauge boson self couplings and a ghost loop, give the well known results for the 
longitudinal and transverse parts of the self energies |||S| 

h , k + k\ 



Ut(K) = -3 6 ab m 2 g (l-^m 
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n T (K) — ^ rn g — 



ko — k J ' 
k 2 \ k k + k 



1 k 2 2k ln h - k 



(4) 

with the thermal SU(2) gauge boson mass m 2 = 2g 2 T 2 /9 and the SU(2) indices a, b of the 
adjoint representation. The remaining two diagrams involving a hard scalar loop show the 
same momentum and energy dependence as in (^) with a contribution to the thermal mass 
of g 2 T 2 /18. In total the expression is changed only by replacing m 2 by tuq = 5g 2 T 2 /18. 
The effective gauge boson propagator in Coulomb gauge is given by 



\ / uu 



k 2 -U L {KY 



Of(K) ^ \ (S„ - f-) Df = ¥ ^ M . (5) 

The Higgs damping rate at zero momentum is defined by 

lH ( p = 0) = ImE* aa (u = m H ,p = 0). (6) 

To leading order in the coupling constant g the imaginary part of the Higgs self energy 
S* comes from the diagram in Fig. 3. The damping mechanism can be read off from this 
diagram by cutting the internal lines [IIJ. The imaginary part of this diagram results 



from the discontinuous part of the effective gauge boson propagator (Landau damping) 
corresponding to the scattering of the soft Higgs off a thermal gauge boson or Higgs via 
the exchange of a soft gauge boson. (There is no pole-pole contribution due to kinematical 
reasons.) Owing to the absence of effective vertices and an imaginary part in the effective 
Higgs propagator there are no bremsstrahlung contributions, i.e. 2 «-> 3 processes, to the 
damping rate in leading order in contrast to the much more complicated damping rate of a 
soft gluon [[J. 

Since we consider the damping rate only at zero momentum, the transverse (magnetic) 
part of the gauge boson propagator does not contribute, leading to an infrared finite result 
due to Debye screening in the effective longitudinal (electric) propagator. 

Using the Matsubara formalism we find 

lH(0) = —-g 2 T dk — — 2— pL(u = u k -m H ,k), (7) 
62n Jo luj: 

where uj\ = k 2 + m 2 H and 

p L {u,k) = --ImD* L {u,k) (8) 

7T 



is the discontinuous part of the longitudinal gauge boson spectral function [ |TTp . 

The remaining integral can be solved only numerically leading to Fig.4, where the Higgs 
damping rate ^njg^T is shown as a function of X/g 2 . The damping rate depends only weakly 



on the Higgs self interaction coupling A, increasing from 0.018 g 2 T at A = to 0.029 g 2 T if 
A tends to infinity. 

The damping rate of the gauge boson at zero momentum is identical to the one calculated 
by Braaten and Pisarski ||. There is no contribution from the polarization diagram by a 
Higgs pair due to kinematical reasons. Furthermore the final result is independent of the 
thermal gauge boson mass mg. (Hence the gluon damping rate does not depend on the 
number of quark flavors in the quark-gluon plasma.) Therefore the gauge boson damping 
rate is given by the result in Ref. || with N c = 2, i.e. 

7l ,t(0) =0.176 g 2 T. (9) 

(The damping rates of a longitudinal and a transverse gauge boson are identical at zero 
momentum.) It should be noted that the gauge boson damping rate is about an order of 
magnitude larger than the one of a Higgs particle. This is caused partly by the fact that 
there is no bremsstrahlung contribution in the Higgs damping and partly by different SU(2) 
group factors. 

Finally we would like to consider the damping rate of a hard Higgs particle with a 
momentum of the order of the temperature or larger. Although it is not related to the 
classical Lyapunov exponent, it might be of interest for cosmo logical problems, since it 
determines the behaviour of thermal and energetic Higgs particles, e.g. relaxation times and 
energy loss, in an electroweak plasma in or near equilibrium. 

The hard damping rate follows from the same diagram (Fig. 3) as the soft one. However, 
now it is sufficient to consider a bare Higgs propagator, i.e. setting tuh = 0. Analogously to 
the case of a hard quark damping rate 1121] we obtain an logarithmically infrared divergent 



result due to the absence of screening in the transverse part of the effective gauge boson prop- 
agator. Assuming an infrared cutoff of the order g 2 T, as for example a magnetic screening 
mass fl3| or by the Bloch-Nordsieck mechanism [I3J we find to logarithmic accuracy 



M P ^T) = ^\n-. (10) 

107T g 
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III. LATTICE SIMULATION 



Now we turn to the computation of the Lyapunov exponent of the coupled Yang-Mills 
Higgs system starting from a classical real-time lattice description. The Hamiltonian of the 
continuum Higgs model, 

h = \eiei + + $ f $ + (A$) f (A$) - ^ 2 $ f $ + A ($ f $) 2 , (11) 



includes the complex Higgs doublet, 



V 



/ 



(12) 



(13) 



-$i + i$ 2 

which can be represented by four real components also as a quaternion 

$ = ($ + *r c $ c ) 
v2 

with c = 1, 2, 3 and the r-s being the Pauli-matrices. These representations are related via 
the scalar product of SU(2) group elements (A, B) = |tr(AB*) 

The Hamiltonian can be approximated on a spatial lattice with lattice spacing a by using 
the variables 



<Px = f *(*) 

as quaternion variables on each lattice site (x = 0, . . . , N 3 — 1) and 

u x,i — c 

as unit length quaternions (([/, [/) = 1) on each lattice link from x to x + i, (i 
The canonically conjugate momenta of these lattice variables are 



(15) 



P =^U ■ 

g 2 



and 



Wx = -o0x, 
^2 



(16) 



0,1,2). 



(17) 



where the dot signals time derivation. The lattice Hamiltonian is [15 



4- ( E(3 - ^) to.. 4>.) + f to.. + E 1 - \v x ,i + (18) 

9 ^ \ x x,i J 

Here mj = a 2 /z 2 and = 4A/g 2 with a being the lattice spacing. We used the following 
lattice approximation of the covariant derivative of the Higgs field 

a ■ (A0x) = U Xt i<f) x+i - 4> x , (19) 

and the identity 

J2(u x ,i<j) x+ i, u X)i (j) x+i ) = Eto*. <f> x ). (20) 

Now scaling the time to if — t/a and the energy to H' = ag 2 H/4 leads to 

h' = e **) + ( 3 - ^) to.. + f to*. 

+ £ i<p x>i , + l - (t^, -T^ + (21) 



X.l 



Here and above the complement link variable, V x> i, is constructed by adding triple products 
of group elements, U, on oriented links closing elementary plaquettes with the chosen link 
(x, i) such that 

E^ tr ^ = E(fW> = \Y,( U *,i> V *J- ( 22 ) 

x,ij x,ij x,i 

Here U Xs ij, the plaquette variable is the product of the four link variables U circumventing 
the plaquette with corner x and lying in the ij plane. 
The scaled equations of motion are 

Ux } i Px,i (frx Ipx 

P x ,i = V Xji + 4) x 4> x+i ip x = W x - F x (23) 



s 



with 



W X = J2 U x ^ x+i + U\^ x -i (24) 

i 

and 

F x = (6 - m f + g f ((j) x ,(j) x )) <p x . (25) 

A gauge transformation which leaves the lattice Hamiltonian (pf) invariant and changes the 
lattice derivative ([19D covariantly is given by 

UL,i = 9 X Ux,i 9x+i 

& = 9l4>* (26) 

with g x being an arbitrary SU(2) group element. The static Noether charge corresponding 
to such time independent gauge transformations becomes 

r * = ^ = E {ul-i,iP x -i,i - P,,iUli) - Ml ( 27 ) 

i 

In order to relate the numerical simulation to physical parameters we consider the average 
energy density 

- 1 V=^&'V (28) 



a 3 N 3 ag 2 a 4 g 2 V3iV 3 

The expression in the last brackets is the scaled energy per link on an iV x N x iV periodic 
lattice. Comparing it with the energy density of an ideal gas of SU(2) gluons and Higgs 
bosons at a fixed temperature T, i.e. with e = 7r 2 T 4 /3 we arrive at 

= y«>T) 4 . (29) 

When obtaining the equivalent temperature the Higgs fields are also considered, since - as 
the results of the numerical simulations - the energy is equally partitioned between electric 
magnetic and rest degrees of freedom in the final, chaotized state. a w = g 2 is the weak 
fine structure constant which can be renormalized at high temperature perturbatively. Its 
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known dependence on the physical temperature (via the perturbative beta function) sets 
the physical scale in the originally scaleless classical simulation. 

Finally we note that in order to ensure the unitarity of the group elements (U, U) = 1 
and its conservation (U,P) = during the solution of the equations of motion Lagrange 
multipliers must be used. Including a small scaled timestep factor dt/a = h into the con- 
jugate momenta and its square into the respective forces we arrive at the following implicit 
recursion scheme, which conserves the static Noether charge |16| 



U' = U + (P' - eU) 

P , = P+(y-fiU + eP') 

<p' = <p + if/ 

if/ = ip + (W — F) (30) 



with fi = (P', P') and e = (U, P'). 



IV. LYAPUNOV EXPONENTS 



We use the following gauge invariant definition of the distance between two different 
gluon and Higgs field configurations, respectively 

D[U,U*]= J £\(U,V) - (U*,V*)\, 

x,i 

£>[$,$*] = J2 !($,$)- ($*,$*)! . (31) 

x,i 

We sampled parallel runs on an N = 10 cubic lattice of randomly initialized U, $ configu- 
rations with unit length. In the dynamical simulation the length of the Higgs field could 
change freely according to the dynamics while the length of the gauge field remained equal 
to one. We repeated each calculation with slightly rotated initial configuration £/*,$*. The 
time evolution of these initially small distances has been observed over a time long enough 
for the saturation. (Distances in a compact space tend to saturate.) 
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As can be seen in Fig. 5 the configurations diverge exponentially in the U as well as in the 
$ configuration space. In the linear range of this evolution the leading Lyapunov exponent 
can be extracted simply 

h=j\nD{t). (32) 

Choosing maximally randomized purely magnetic gauge field configurations and similar 
Higgs field configurations initially we studied the chaotic behavior due to the extraction of 
leading Lyapunov exponents for different values of gj = 4X/g 2 . As a result of this study 
the Higgs and gauge boson Lyapunov exponents are plotted in Fig. 6. As expected, the 
Lyapunov exponent of the gauge boson does not depend on A and agrees with twice of the 
damping rate (||) within the numerical accuracy. 

Surprisingly the Higgs Lyapunov Exponent shows a behavior qualitatively and quantita- 
tively different from that of the damping rate obtained in the resummed one loop approxi- 
mation (Fig.4). We remind, that this is contrary to the case of the pure gauge field system, 
where the Lyapunov exponent has been found to be twice of the gluon damping rate at zero 
momentum ||. 

The result of the classical lattice dynamics can be understood easily: for very high values 
of A the $-fields are constrained to keep a rather definite value of ($,$). The number of 
light degrees of freedom is reduced from four to three, so the Higgs sector becomes less 
chaotic. In the simulations presented here the gauge fields were initialized with random 
phases covering the whole group SU(2) uniformly in the magnetic sector, so these degrees of 
freedom could not take any energy from the Higgs fields over. The electric sector was able to 
receive energy also from the Higgs, but because of the final equipartition with the magnetic 
degrees of freedom the total energy converted into chaotic motion remained limited. 

From the results of the numerical simulation presented here we must conclude that the 
Higgs damping rate, contrary to the gluon damping, is not closely related to the chaotic 
behavior of the classical fields. While the chaotic behavior of Higgs fields is rather set by the 
damping rate of gluons which couple to the Higgs covariantly and decreases somewhat at 
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high A values, their resummed damping rate at high temperature increases with increasing 
Higgs-self coupling due to the thermal Higgs-mass which sets the scale in that calculation 
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FIGURES 




FIG. 1. HTL contributions to the Higgs self energy. 




FIG. 3. Higgs self energy contribution to the damping rate. 
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FIG. 4. Higgs damping rate 7#/<7 2 T at zero momentum as a function of X/g 2 ■ 
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FIG. 5. Time evolution of the distance of two initially adjacent Higgs field configurations at 
A\/g 2 = 1024. 
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FIG. 6. The gauge field (upper curve) and Higgs field (lower curve) leading Lyapunov exponents 
as a function of 4A/g 2 at high temperature. The numerical results have an error bar of about 20%. 
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